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Abstract 

. This article is the first in a series deahng with the thermodynamic 

' properties of quantum Coulomb systems. 

. In this first part, we consider a general real- valued function E de- 

fined on all bounded open sets of R'^. Our aim is to give sufficient 
conditions such that E has a thermodynamic limit. This means that 
' the limit E{Vtn)\Q,n\~^ exists for all 'regular enough' sequence 0„ with 

growing volume, |r2„| oo, and is independent of the considered se- 
quence. 

The sufficient conditions presented in our work all have a clear phys- 
ical interpretation. In the next paper, we show that the free energies of 
many different quantum Coulomb systems satisfy these assumptions, 
hence have a thermodynamic limit. 



1 



Contents 



Introduction 2 

1 General Theory of Thermodynamic Limit 9 

1.1 Definitions 9 

1.2 Thermodynamic limit for a special sequence of sets 11 

1.3 Thermodynamic limit for general sets 15 

2 Proofs 21 

2.1 Proof of Theorem 1: limit for the reference set 21 

2.2 Proof of Theorem 2: limit for general domains 25 

2.3 Proof of Proposition 2: regularity of inner approximation . . 36 

2.4 Proof of Proposition 3: (A6) implies (A5) 37 

2.5 Proof of Lemma 1: regularity of open convex sets 38 

References 38 



Introduction 

Ordinary matter is composed of electrons (negatively charged) and nuclei 
(positively charged) interacting via Coulomb forces. The potential between 
two such particles of charges z and z' located at x and x' in is 



\x — x'\ 

There are two difficulties which occur when trying to describe systems com- 
posed of electrons and nuclei. Both have to do with the physical problem of 
stability of quantum systems. 

The first is due to the singularity of l/|x| at 0: it is necessary to explain 
why a particle will not rush to a particle of the opposite charge. One of the 
first major triumphs of the theory of quantum mechanics is the explanation 
it gives of the stability of the hydrogen atom (and the complete description 
of its spectrum) and of other microscopic quantum Coulomb systems. 

Quantum Mechanics or more precisely the uncertainty principle explains 
not only the stability of tiny microscopic objects, but also the stability 
of gigantic stellar objects such as white dwarfs. Chandrasekhar's famous 
theory on the stability of white dwarfs required, however, not only the usual 
uncertainty principle, but also the Pauli exclusion principle for the fermionic 
electrons. 
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Both the stabihty of hydrogen and the stabihty of white dwarfs simply 
mean that the total energy of the system cannot be arbitrarily negative. 
If there was no such lower bound to the energy one would have a system 
from which it would be possible in principle to extract an infinite amount of 
energy. One often refers to this kind of stability as stability of the first kind 
[15, 16]. If we denote by E{N) the ground state energy of the system under 
consideration, for N particles, stability of the first kind can be written 



In proving (1) for Coulomb systems, a major role is played by the uncer- 
tainty principle which is mathematically expressed by the critical Sobolev 
embedding H^{M.^) ^ L^(R^). The latter allows to prove Kato's inequality 



for all e > 0, which means that the Coulomb potential is controlled by the 
kinetic energy. 

The second issue concerns the slow decay of l/\x\ at infinity and this 
has to do with the macroscopic behavior of quantum Coulomb systems. It 
is indeed necessary to explain how a very large number of electrons and 
nuclei can stay bounded together to form macroscopic systems, although 
each particle interacts with a lot of other charged particles due to the long 
tail of the Coulomb interaction potential. Whereas both the stability of 
atoms and the stability of white dwarfs were early triumphs of quantum 
mechanics it, surprisingly, took nearly forty years before the question of 
stability of everyday macroscopic objects was even raised (see Fisher and 
Ruelle [9]). The rigorous answer to the question came shortly thereafter in 
what came to be known as the Theorem on Stability of Matter proved first 
by Dyson and Lenard [5, 6]. 

The main question is how the lowest possible energy E{N) appearing in 
(1) depends on the (macroscopic) number N of particles in the object. More 
precisely, one is interested in proving a behavior of the form 



This behavior as the number of particles grows is mandatory to explain why 
matter does not collapse or explode in the thermodynamic limit. Assume 

that (2) does not hold and that for instance E{N) ^n^oo cN^ with p ^ 
1. Then \E{2N) — 2E[N)\ becomes very large as N ^ 1. Depending on 



E{N) > -GO. 



(1) 



(2) 
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p and the sign of the constant c, a very large amount of energy will be 
either released when two identical systems are put together, or necessary 
to assemble them. The constant e in (2) is interpreted as the energy per 
particle. 

Stability of Matter is itself a necessary first step towards a proof of (2) 
as it can be expressed by the lower bound: 

E{N) > -kN. (3) 

Put differently, the lowest possible energy calculated per particle cannot be 
arbitrarily negative as the number of particles increases. This is also often 
referred to as stability of the second kind [15, 16]. 

A maybe more intuitive notion of stability would be to ask for the volume 
occupied by a macroscopic object. Or more precisely, what is the volume 
of the object when its total energy is close to the lowest possible. It can be 
shown in general that this properly defined volume is proportional to the 
number of particles A^. Denoting by a domain in which is occupied 
by the system under consideration and by E('n) its (lowest possible) energy, 
(2) then reads 

m) -mHoo e\n\ (4) 

where \ is the volume of fl. Similarly, stability of the second kind is then 
expressed as 

E{n) > -K\n\. (5) 

Large quantum Coulomb systems have been the object of an important 
investigation in the last decades and many techniques have been developed. 
A result like (3) (or equivalently (5)) was first proved for quantum electrons 
and nuclei by Dyson and Lenard [5, 6]. After the original proof by Dyson and 
Lenard several other proofs were given. Lieb and Thirring [22] in particular 
presented an elegant and simple proof relying on an uncertainty principle 
for fermions. The different techniques and results concerning stability of 
matter were reviewed for instance in [15, 16, 17, 23, 27]. 

It is very important that the negatively charged particles (the electrons) 
are fermions. It was discovered by Dyson in [4] that the Pauli exclusion prin- 
ciple is essential for stability of Coulomb systems: charged bosons are alone 
not stable because their ground state energy satisfies E{N) ~ -CN^^, see 
[4, 2, 21, 26]. 

A result like (2) (or equivalently (4)) was first proved for quantum 
Coulomb systems by Lieb and Lebowitz in [18] for a system containing 
electrons and nuclei both considered as dynamic quantum particles, hence 
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in particular invariant under rotations. Later Fefferman gave a different 
proof [7] for the case where the nuclei arc fixed particles placed on a lattice, 
a system which is not invariant under rotations. 

Instead of the ground state energy, one can similarly consider the free 
energy F{Q,P,fi) at temperature T = 1/(3 and chemical potential ji. A 
precise definition of this quantity will be provided for the three examples 
that we shall consider in [13]. One is interested in proving the equivalent of 
(4) 

F{n,l3,ij)^\a\^^f{f3,lj)M (6) 

where fi) is the free energy per unit volume. The pressure is then given 
byp(/3,/x) = -/3/(/3,/i). 

In this work, we provide a new insight in the study of the thermodynamic 
limit of quantum systems, by giving a general proof of (4) or (6) which can 
be applied to many different quantum systems including those studied by 
Lieb and Lcbowitz in [18] or Fefferman in [7]. Our goal was to identify the 
main physical properties of the free energy which are sufficient to prove the 
existence of the thermodynamic limit. For this reason, this first paper will 
be dedicated to the study of an abstract energy 

VL ^ E{n) 

defined on all open bounded subsets of with values in M, for which the 
limit 

lim -TT^ = e (7) 

n— >oo \iln\ 

holds as \Qn\ oo. In principle, E can represent the ground state energy 
or the free energy at temperature T > of any quantum (or even classical) 
system. 

We shall apply our general framework to several quantum Coulomb sys- 
tems in another paper^ [13]. More precisely, we shall treat in [13] three 
different systems: 

• the quantum crystal for which the nuclei are fixed particles on a lattice, 
with a possible local rearrangement and local defects; 

• dynamic quantum nuclei and electrons in a periodic magnetic field; 

• optimized classical nuclei with quantum electrons. 



^Equations or results with reference n in the second paper [13] will be denoted as Il.n. 
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Like in previous works, our method consists in first showing the existence 
of the Hmit (7) for a specific domain A which is dilated (and possibly rotated 
and translated). Usually A is chosen to be a ball, a cube or a tetrahedron. 
In the applications [13] we always use a tetrahedron as we shall heavily rely 
on an inequality proved by Graf and Schenker [11] for this type of domains, 
see Proposition II. 2. 

The second step consists in showing the existence of the limit (7) for 
any (reasonable) sequence of domains {O^} such that jQ^j — > oo. This is 
important as in principle the limit could depend on the chosen sequence, a 
fact that we want to exclude for our systems. Here we do not specify what 
a "reasonable" sequence is but some properties will need to be assumed to 
ensure that boundary effects always stay negligible. It is not very surprising 
that proving the existence of the limit for any sequence {fin} is much more 
complicated (hence requires more assumptions on the energy E) than for 
a simple sequence made of the reference set A. For this reason, we have 
divided our results in two parts and first treat the simpler case of A. 

It is to be noticed that our method (relying on the Graf-Schenker in- 
equality) is primarily devoted to the study of quantum systems interacting 
through Coulomb forces. It might be applicable to other interactions but 
we shall not address this question here. 

We now describe vaguely the assumptions which we found to be sufficient 
on the function E to prove the existence of the limit (7). Of course we shall 
give more details later on. 

To prove the existence of the limit for sequences made of the reference 
set A, we need five main assumptions. In principle, A can be any convex 
bounded open set. The first is the normalization condition: 

(Al) E{^) = 0. 

The second important assumption on E is the stability of matter: 

(A2) vn G M, E{n) > -K\n\. 

as we explained before. In the following Ai is by definition the set of all 
bounded open subsets of whereas TZ will be a subclass of 'regular' sets in 
a sense which will be made precise later. Determining the correct regularity 
assumption is indeed a subtle task which we do not explain more in this 
introduction. We also need some translation invariance property: we assume 
that the following limit exists for all ^ & TZ: 
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Examples are given by fully-translation invariant systems, or periodic sys- 
tems. But we can also treat perturbations of those. The following continuity 
property will play an important role: 

(A4) yn', nen with n' c n, E{n) < E{n') + K\n \ n'\ + error. 

It essentially says that a small decrease of will not decrease too much 
the energy. A similar property was used and proved in the crystal case by 
Fefferman [7, Lemma 2]. In the above formula (and below), 'error' describes 
an error term which will be precised later. Indeed, as we will see the specific 
form of the error terms will play an important role in the proof of the 
existence of the thermodynamic limit. Finally, our method is mainly based 
on the following inequality 



which compares the energy of 0, with the energy of the reference set A, 
averaged over rotations and translations of A inside 17, see Figure 1. This 
type of property was first remarked and used by Conlon, Lieb and Yau [2, 3], 
for quantum systems interacting with the Yukawa potential and A being a 
cube. For Coulomb interactions, it was proved by Graf and Schenker [11, 10] 
in which case A is chosen to be a tetrahedron. Fefferman also used a proof 
by means of a lower bound in his study of the crystal case in [7]. This will 
be explained in details in [13]. 



Figure 1: Assumption (A5): the energy of 0, is greater than the energy per 
unit volume of the reference set A, averaged over all possible translations 
and rotations inside il. 

It will be shown in Theorem 1 that assumptions (Al)— (A5) are suffi- 
cient to obtain the existence of the limit 





lim 



oo 



= e 
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for all sequences in — > oo, {{Rn,Un)} C 50(3) x M^. 

In order to prove the existence of the limit for any sequence we 
need the additional assumption that A is a polyhedron which can be used 
to build a tiling of M'^ and the following new condition. For all O E 7^, we 
assume that there exist numbers /j^'"^ and a function 5 : 7W ^ M such that 
the energy can be decomposed in the form^ 

( A6) E{n) = Y,e( nn Af + i 4""'^^ " ^(^) 

+ ^SmnA5^'''M + error 

where {A^^'"^} are all the elements of the tiling, rotated by ii G 50(3) and 
translated by u G , see Figure 2. The number /-^'"^ is interpreted as the 

interaction energy between the sets A^^'"^ and a[.^'"'*. The numbers Ij^'^^ 
must satisfy some properties which we do not mention in this introduction, 
for this interpretation to be meaningful. 




Figure 2: Assumption (A6): the energy in the big set Q can be decomposed 
into energies in small sets A^^'"^ of the tiling plus the interaction between 
them and the difference between entropies. 

The assumption which corresponds to (A5) is that the total interaction 
energy is nonnegative when averaged over translations and rotations of the 
tiling (the numbering of the equations corresponds to that of Section 1.3): 




—error. 



"^In the text we indeed introduce two assumptions (A6.1) and (A6.2) corresponding 
to, respectively, a lower and an upper bound on the error term. 
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Finally, the function S{Cl) appearing in (A6) is assumed to satisfy certain 
properties similar to the ones of entropies in statistical mechanics. In the 
applications [13], S will indeed just be the usual quantum entropy. The 
main assumptions'^ on S are its normalization 

(A6.5) 5(0) = 

and the strong subadditivity property: 

(A6.6) s{ni u O2 u Qs) + 5(02) < s{ni u Q2) + s{n2 u O3) 

for all disjoint subsets Oi, ^2 and ^3. 

Conjectured by Lanford and Robinson [14] the strong subadditivity (SSA) 
of the quantum mechanical entropy was proved by Lieb and Ruskai in 
[19, 20]. Robinson and Ruelle noticed in [24] (see also Wehrl [28]) that 
SSA was a necessary property in order to define the entropy per unit vol- 
ume of translation-invariant states: they used it to control boundary terms 
of any arbitrary Van Hove sequence of domains. To our knowledge, SSA 
was never used in the way we do it to prove the thermodynamic limit for 
interacting systems. 

The paper is organized as follows. In Section 1 we state our main results, 
starting from the case of a special sequence Qn = Rn^n^ + "Un for which 
we need less assumptions than for the general case. Proofs are given in 
Section 2. Let us mention that the results of this paper and of [13] have 
been summarized in [12]. 

Acknowledgment. M.L. acknowledges support from the ANR project 
"ACCQUAREL" of the French ministry of research. 

1 General Theory of Thermodynamic Limit 

1.1 Definitions 

Regularity property of sets in 

Let us first define the class of sets for which we shall be able to prove the 
thermodynamic limit. In the whole paper, we denote by \A\ the Lebesgue 

^In the text only the diflference S'(fi) — SiQ, fl A^^'"') is considered. See the details 
in Section 1.3. 
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measure of any measurable set ^ C R^, and by #B the cardinaUty of any 
finite set B. Let 

is open and bounded} . 

We shall restrict ourselves to the three dimensional case for simplicity but all 
the results of this first paper can easily be generalized to the N dimensional 
case. We shall need the following 

Definition 1 (Sets with regular boundary). Let r] : [0,to) i-^ [0, oo) be a 
real function^ with limt_»o r]{t) = 0. We say that Q E M has an ry-regular 
boundary in the sense of Fisher if 

Vt G [0,to), {x G I d{x,dn) < < Mv{t), 

where d^} = f7 \ $7 is the boundary of ft . 

We denote by 7^^ C M. the set of all subsets Vt ^ M. having an rj-regular 
boundary in the sense of Fisher. 

The definition of r/-regular boundary was introduced by Fisher [8, page 
394] (see also [18, page 351]) and it is invariant by scaling. In the following, 
we shall not consider any function 77, but for simplicity only those belonging 
to the following class 

:= |r? : [0,c) [0,oo) | r}{t) = at^, a > 0, c> 0, be (0,1]| . 

Most of our results could be easily generalized to a more general class of 
functions rj. But the class of sets having Fisher's boundary property for 
some ?7 G f already contains very "perturbed" sets for which the area of the 
boundary can be of order greater than irip^^, see the footnote of [8, page 
394]. 

Notice the following simple result, whose proof is given in Section 2.5. 

Lemma 1. Any open and bounded convex subset of belongs to TZrt for 
all r]{t) = at, t & [0,^0)7 ''JJith a large enough and to small enough. 



The sliding group 

We consider the group consisting of orientation preserving isometrics of M^, 
i.e. the semidiroct product G = X SO{3). This subgroup G of the 

''in the whole paper, we use the convention that 77 is given with its domain of definition, 
i.e. with to > 0. 
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Euclidian group (it is sometimes called the Special Euclidian group) acts on 
as follows: {u, R) ■ x = Rx + u. Let dX denote the Haar measure of G. 
Then 

/ f{gy)dX{g) = [ f{x)dx (8) 
Jg Jk.^ 

for any / G L^(R^) and any y G M^. 

Of course, G also acts on M and it stabilizes 7?.^ C M for any rj e S. 

1.2 Thermodynamic limit for a special sequence of sets 

Let us consider a function E : M ^ W. Our main goal is to prove the 
existence of a thermodynamic limit under suitable assumptions on E. This 
means that there exists a real number e such that 

Ejnn) _ 

iim = e 

for any 'regular' sequence of sets {^n}n>i with |r2„| — > oo as n — oo. 

The regularity of sets will be measured differently depending on the 
applications. For the general theory, let us consider some fixed function 
r] E S defined on [0, c), and an abstract class of sets TZ Q M.. Here we only 
assume that Tl is invariant under translations, rotations and scaling, and 
that it contains sets having at least Fisher's 77-regularity boundary property: 

(PI) One has 7^ C 7^^. If J2 g 7^, then gm G 7^ for any £ > 1 and any 
geG. Moreover, G 7^. 

For the crystal case, Ti will consist of the sets having the cone property and 
a regular boundary, see [13]. In some other cases, we shall be able to take 
TZ = TZff for some fixed rj E S. We also define the regularized volume of any 
Q e M as follows 

\n\r := mi{\ni nen, non}, (9) 

i.e. \Q.\r is the smallest volume of all regular sets Cl eTZ containing $7. Notice 
the obvious inequality > 

Also we shall need another class TZ' of subsets of M^. We only assume 
at this stage that there exists some rj' G £ defined on [0, c') with rj < rj' on 
[0, min(c, c')) such that 

(P2) 7^ c 7^' c 7^^,. 
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The role of Tl' will appear more clearly in the next section. For this section, 
the reader can simply take TZ' = TZ. 

In this section, we start with the problem to prove the existence of the 
thermodynamic limit for specific sequences made of a reference set A £ TZ, 
i.e. J7„ = gn^n^ with {gn\ C G and — > oo. We require adequate 
properties on the function which will need to be strengthened later to 
obtain the limit for general sequences. 

For the moment, we only assume that the reference set A satisfies the 
following property: 

(P3) A G 7^ is a convex set such that € A. 

Thus (1 — A) A C A for any < A < 1. Later, we shall need to assume 
that A can be used to define a tiling of M^, but this is not necessary at this 
point. 

For the function we assume that there exists a subset C TW, a 
function e : 7?, ^ M, two positive constants 5, k > and a function a with 
lim^^oo ct(^) = such that the following hold for ^>l•. 

(Al) (Normalization). E{^) = 0. 

(A2) (Stability). For all Q e M, we have 

E{n) > -K\n\. 

(A3) (Translation Invariance in Average). For all ^ E TZ, 



|S(0,L)|-i / E^pldu-e{n) 

Jb(o.l) m 



< a{L), (10) 



'B(0,L) 

where B{0, L) denotes the ball of center and radius L. 

(A4) (Continuity). For all G 7^ and Vt' G 7^' with C O and d(5J^, dVt') > 
6, we have 

E{n) < E{n') + K\n \ n'\ + \n\ai\n\). 

(A5) (Subaverage Property). For all Q G M5, we have 

1 _ Q,(l) f 

m) > la^^^^^' ^^^^^ '^^^^^ ~ '^'"^ ^^^^ 
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Remark 1. Notice that (Al), (A2) and (A4) imply that there exists a 
constant k' > such that 



(12) 



for all € 7^ and for large enough (take $7' = in (A4)). However the 
upper bound need not be truc^ for all Q, in Ai. 

Remark 2. Assumption (A3) is a generalization of the periodicity require- 
ment. Examples of functions E satisfying (A3) contain 

1. functions satisfying (12) and which are A-periodic for some discrete 
subgroup A C R-^ with compact fundamental domain R^/A, hence in 
particular fully translation-invariant systems, V-u G M^, E{Q, + u) = 



2. finite sums of periodic functions as above; 

3. functions which are almost A-periodic, uniformly with respect to Q. 

In the proof of our results, we actually do not need (A3) for all G 7^, but 
only for Q = giA with g & G and i large enough. 

Remark 3. The most important assumption is of course (A5) which allows 
to bound from below the energy in a large set fl by an average of energies 
of the (smaller) scaled reference set £A, as explained in introduction. We 
introduced the set for which property (A5) is true. In this section we 
shall simply take = Ma with 



for some fixed ^ 1; but in the limit for general sets in the next section, we 
shall need (A5) for all domains J7 (i.e. we will choose M5 = M). Notice in 
all cases the error is assumed to depend on \^\r, that is to say on the volume 
of the smallest regular set containing fl. 

Remark 4. Property (A4) is also very important. It essentially says that 
a small decrease of Q will not decrease too much the energy. A similar 
property was used and proved in the crystal case by Fefferman [7, Lemma 
2]. In the proof of the existence of the thermodynamic limit for general 

®In the crystal case studied in [13], the upper bound does not hold for all f2 € A^. For 
the other two examples of [13] the upper bound is true, more precisely one has E{Q,) < 
for any Q in M. 



E{ny, 



Ma ■■= {gi^, geG, £>io} 



(13) 
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domains, we shall need this property for domains $7' which are a bit less 
regular than sets in TZ (see Section 1.3), explaining the introduction of the 
class TZ'. For the result of this section, we could simply take TZ! = TZ but 
later on we shall need to cover the case TZ' 2 ^• 

Remark 5. Notice that the function g ^ E(^D, D g ■ (^A)) has compact 
support by (Al) and is bounded below by (A2): 

E{nng-{eA)) > -K\n\. 

Therefore (A5) in particular implies that g E(p,r] g ■ (^A)) belongs to 
L^{G,dX) for any U e M5. 

Remark 6. Assume that is a local functional of the form 

E{n)= f x{x)dx (14) 
Jn 

for some fixed function x £ (I^^ j ^) satisfying 

\BiO,L)\'^ [ xix + u)du-Xoc{x) 
Jb(o.l) 



lim 

L—*oo 



= (15) 



for some function Xoo : M 1— > M, uniformly with respect to x (it can be proved 
that Xoo is necessarily constant, see Lemma 4 below). See [1] for comments 
on the property (15). Then (A1)-(A5) are satisfied for any 0, £ A4 and 
any reference set A. In particular (A5) is always an equality, as 

[ E{nng-A)dX{g) = ^ [ dx [ dX{g)ln{x)lA{9-^x)x{x) = E{n) 
\^\ Jg \^\ JR3 Jg 

(16) 

for any set A e M, and where we have used that 

dX{g)lA{g-^x) = \A\, 



IG 

see (8). For the volume E{p) = x = 1, this in particular gives 



If a function E takes the form (14) for some x £ L°°(M , M) satisfying (15), 
then the thermodynamic limit is easily proved to be the constant Xoo- 
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We are now able to state our 



Theorem 1 (Thermodynamic limit for the reference set A). We assume 
that TZ, TZ' and A satisfy (P1)-(P3) and that the function E satisfies 
(A1)-(A5) with = A^a defined in (13). We denote 

= -JET 

forgeG = R^x SO{3) andl>l. 

Then there exists a real constant e such that e^ converges uniformly to 
e on G as i ^ oo. In particular, for all sequences {£n}n>i ^ [l^cc) with 
in ^ oo and {gn}n>i Q G, 

iim — = e. 

n— *oo |-t^A| 

Additionally, the limit e does not depend on the reference set A ; if the 
above properties hold true for another set A' and if = A^a U A^a'; then 
E(^g£A'^\£A'\~^ converges to the same limit e as £ ^ 00, uniformly on G. 

The proof of Theorem 1 is given in Section 2.1 below. 

Remark 7. We notice that Theorem 1 holds true if we replace (A2) by the 

weaker assumption involving the regularized volume 

(A2') (Weak Stability). For any Q e M, then 

E{n) > -K\n\r. 

However (A2) will be needed in the next section where we tackle the case 
of a general sequence of domains. 



1.3 Thermodynamic limit for general sets 

In Theorem 1, we have stated the existence of the thermodynamic limit for 
domains A satisfying certain assumptions. However in the applications it is 
not clear for which domain A, property (A5) will hold true. At present, for 
Coulomb systems an inequality of the form ( A5) is only known for simplices 
[11, 13]. Hence, in order to prove the existence of the thermodynamic limit 
for any (regular) sequence of sets, we need to add some assumptions. The 
idea will be to decompose any domain as a union of (translated and rotated) 
reference sets A. For this reason we will now assume that A can be used 
to build a tiling of the whole space, and we will require some localization 
properties on the energy E. First, we introduce the following 
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Definition 2 (Tilings of R^). LetT be a discrete subgroup ofG with compact 
quotient G/F. We say that the open set A C defines a F-tiling of if 

• /xA n ii'A = $ if fj,,ii' eV and jJL ^ ji' ; 

In the following, we assume additionally to (P1)-(P3) that 

(P4) A is a polyhedron which defines a F-tiling for some discrete subgroup 
r of G with G/T compact. 

In the Coulomb case [11, 10], we shall essentially (up to a translation) take 
r = X O, where O is the group of order 24 consisting of the pure rotations 
of the octahedral (or hexahedral) group, i.e. the symmetry group of the 
cube. The set A is the simplex which makes up a 24th of the cube. 

Remark 8. Notice that there is no loss of generality in assuming that the 
convex set A is a polyhedron: it can easily be proved by means of the finite- 
dimensional Hahn-Banach Theorem that any convex set defining a T-tiling 
is necessarily a polyhedron. 

Given A which defines a F-tiling, any regular set Vi ^TZ can be approx- 
imated by unions of elements of the (dilated and rotated) tiling, see Figure 
3. The following proposition tells us that the so-approximated set has a 
boundary which stays regular, although it might be a bit less regular than 
that of Q,. It will be very useful in the proof of our results. 

Proposition 2 (Regularity of the boundary of the inner approximation). 
Let (5, To > 0, f2 € and introduce for any g E G, T,i > 

^T,M = ^5/^(1 + 

If 

ArA9m= U Ar,eMCn, (18) 

AtGr, 
Ar,«,9(/i)cn, 

d{dAr,e,g{ix),dn)>S 

then there exists three constants io > 0, £i > and m> 1 such that 

I^^Ar/,g{n) I Q G 7^r,, g€G, < T < TO, io<i< l^^l^-^^^i} C TZfi, 

with fj{t) = mrjit), t G [0, c/m) where we recall that the domain of definition 
off] is [0,c). 
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Remark 9. Since 6 will be fixed in accordance with (A4), we do not ex- 
plicitely emphasize that AT-/^g{^) depends on 5. 

Proposition 2 is shown in Section 2.3. In the proof of our result, we 
shall need to know that the energy E[^) does not vary too much when J7 is 
replaced by its inner approximation = ^,-^^ ^(0). This will be true pro- 
vided we can apply (A4), leading to the natural assumption that Ar^t^g{^) 
belongs to TZ' for any £ TZ (the constant 6 is chosen like in (A4)): 

(P5) For any tq > 0, there exists two constants > 0, > such that 

\^Ar,i,g{n) \nen, g^G, o < r < to, £o<i< m'^^^h} e 

We need to add some assumptions on the energy E to prove the existence 
of the thermodynamic limit for general sets. We assume that there exists a 
class of domains C Ai such that the following holds 

(A6) (Local Energy Decomposition) . For all 17 G C Ai and all ^ > 1 
there exist maps 

E^ -.G^ R, 
I? ■■ {{9,9') G G X G I g-'g' e F \ {0}} ^ R 
sf :G X {V CT} ^R 
such that the following conditions hold for all g G: 
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(A6.1) (Lower hound). 



E{n)>Y,E^{9t^)+^ E I?{9f^,9^)-sf{9,r)-\n\a{£). (19) 



(A6.2) fC/pper hound). For all P C T, 




- s^^ia.V) + O n U +t(^))A a(£), 



for some continuous function r : (0, oo) — [0, oo) which does not 
depend on f2 and satisfies lim^_»oo r(£) = 0. 

(A6.3) (Interaction Average). 



(A6.4) (Interaction Support). iy{gii,gu) = if egfi{l + T{i))Ann = 
or if £c/i/(l + t(£))A n = 0. 

(A6.5) (Normalization of sf). For any /x € F and g & G, 



(A6.6) (Strong Subadditivity of s^). For any disjoint subsets Vi, V2, 



Let us make some comments. Assumption (A6) essentially means that 
our energy arises from a physical system containing at most two-body inter- 
actions. The reader has to imagine that Ef^gji) represents the total (free) 
energy in each £giJ,A, whereas lf'{gii,gv) describes the interaction energy 
between ig/iA and igi/A. Our main result below can easily be extended to 
the case of a /c-body interaction but we do not give further details as this 
would complicate (A6). 




per 



sfig,{n})=sf{g,ili)=0. 



^3 and any g E G, 



4\g, Vs) + 4\g, P2) < 4\g, Vi U P2) + s\\g, V2 U V3). 
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Assumption (A6.3) means that averaging over translations and rota- 
tions of the tihng yields a global interaction energy which is essentially 
nonnegative. This is a specified aver aging- type property which is more pre- 
cise than our assumption (A5) above. Indeed, we shall prove that (A6) is 
stronger than (A5). This is expressed in the following result, whose proof 
will be given in Section 2.4. 

Proposition 3. Let TZ, TZ' and A satisfy (P1)-(P5), and E satisfying 
(A1)-(A4). If E satisfies (A6) for all Vt G Mq where Mq C M, then E 
also satifies the averaging property (A5) for all € M.%. This means we 
can assume that A^e ^ -^5 

In our study, we shall need (A5) for all domains J7 G whereas we 
only use (A6) for Q. G A^a which was defined in (13). If (A6) is true for 
all domains O G A(, we do not need assumption (A5) by Proposition 3. 

In the applications, the function will be related to the entropy of the 
system. Formally, it will be the difference between the entropy in the big 
domain and the entropies of the small sets: 



For the original tiling {IgjiA] ^j^^y , it might be difficult to prove that the 
energy behaves as stated, in particular because of the considered boundary 

conditions and localization issues. For this reason, one may need that each of 
the sets used for the tiling slightly overlap, which is the role of the correction 
r(^). Notice that since G A, then A C (1 -|- r(^))A. 

Remark 10. For the sake of simplicity, we have used the same function 
a as in the previous section. We recall that lim^^ooa(^) = 0. Symetriz- 
ing the above properties if necessary, we can always assume that if {■■,■) is 
symmetric: 



Note that it would be natural to assume in addition to (A6.4) that E{gii) = 
when £gn{l + T{i)) Ann = and that s^{g,V) = when £5/x(l-|-r(^))An 
O = for any n EP, but we actually do not need it in the proof. 

Remark 11. The stability condition (A2) is used in the proof to control 
some complicated boundary terms. Our proof does not work under the 
weaker assumption (A2') introduced in the previous section. 




yg G G, V/x, 1/ G r, fij^ V, I^ {gji, gv) = I^ {gv, gfi). 
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Remark 12. Taking 7^2 = in (A6.6) and using (A6.5), we deduce that 
sp is subadditive: 

s'lig^Vi U V^) < s2ig,Vi) + s?{g,P3). (20) 

Then, using that sf{g,{iJ,}) = for any G F by (A6.5), we also deduce 
by induction that s is monotone 

Pi C ^2 C r ^ s^ig,r2) < sf{g,Pi), (21) 

hence that 

VPCr, sf{g,r)<0 (22) 

for any g E G. 

Remark 13. In the proof we do not use the fuh strong subadditivity 
(A6.6), but rather the weaker property 

sng,r)<^ ^ s2{g,{nM) (23) 

which is satisfied for any strongly subadditive function satisfying (A6.5), 
as proved in Lemma 9. However the generalization to a fc-body interac- 
tion would require an estimate different from (23) and for this reason we 
prefer to keep the more general assumption (A6.6). Also it is well-known 
that the strong subadditivity of the entropy is a useful tool in the study of 
thermodynamic limits [24, 28]. 

We are now able to state our main 

Theorem 2 (Thermodynamic limit for general domains). We assume that 
TZ, TZ' and A satisfy (P1)-(P5) and that the function E satisfies (Al)- 
(A6) with = M and Mf, = Ma- 

Then, for any sequence {^n}n>i ^ ^ such that diam(ri„)|r2„|~^/^ is 
bounded and — ^ oo, we have 

n— >oo 

where e is the limit obtained in Theorem 1. 

We recall that Ma was defined in (13). The rest of the paper is devoted 
to the proof of our results. Applications are given in the second part [13]. 

Remark 14. As we have TZ C 7?.^, we know from [18, Appendix A p. 385] 
and [8, Lemma 1] that if each set f2„ is connected, then automatically 
diam(0n)|ri„|-i/3 < c for all n. 
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2 Proofs 

2.1 Proof of Theorem 1: limit for the reference set 

Notice that applying (12) to g£A which is in TZ by (PI) and (P3), we infer 

that €£ is uniformly bounded on G. 

Step 1. Translation invariance of e defined in (A3). We start by proving 
the following 

Lemma 4. For allQ. ^TZ and all u G M? , one has e{Q, + u) = e(0). 
Proof. Let O G 7?. and G M^. Then, for L large enough, 



Jb{o,l) I"I Jb{u,l) 



\B{0,L)\Biu,L)\ + \Biu,L)\BiO,L)\ 
- ^ \B{Oj:)\ 

by (12) and since Q + w G 7?. for any G M'^ by (PI). The result is obtained 
by taking the limit L — oo. □ 

Step 2. Lower bound. We now prove the important 

Lemma 5 (Bound from below). There exist a decreasing function ai tend- 
ing to zero at infinity and a constant Ki > 1 such that, for all fj E £ ( defined 
on [0, c') ), all ^ €TZfjn M5 and all i>l, 

^(")>er-Sa.W-<:^(^) (24) 



\n\ - ' ^ 

provided £\Q\~^/^ < c' /ki, and with 



ef := / e{mA)dR. (25) 
JsO(3) 

Proof. Let G T^r] H A^s and i > 1. We introduce the set of 'interior points' 

Ie = {ueM.^ \ {MA + u)^9,, MRe SO{?>)}. 
and of 'boundary points' 

Bt = {ue^^ \ 3Re 50(3), {RIA + u) n O 7^ 0} \ I^. 
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Notice that since G A, then 1^ C $7. By (A5), we obtain 

E{n) > {i-a{e))[ dx{g)^j^ - 

JlexSOis) 1*^1 

<iAto)^(5^«l (20) 

> (1 - a{i)) [ duf du^^^^^±^ - MM^) 

JS0{3) Jli \^^\ 

-{1 - a{e))K\Bi\. (27) 

where we have also used (A2). 

We now use that (A5) is an equality for the local function x '■ u i-^ 
E {RiA + u) \£A\~^ and with the ball 5(0,1) := A as reference set, as 
explained in (16). We obtain 

Jie I^A| \B{o,e)\JG ^^'Ji.ngBm m 

1 f , f E{RIA + a) 

" \mw\L LnBi.,e) m ^ ' 

We introduce similarly to I( and B^, 

/; := {x G I B{x,e) C I^}, 

B^:= {xeR^l B{xJ)r]n^(li}\4. 

Notice that I[ Q h Q ^ and C Ba(. for some constant a depending only 
on A. Then (28) becomes by (A2) 

f ^ E{RlA + u) ^ f ^ 1 f E{RiA + u) ^ 
Using (A3) and Lemma 4, we obtain that 

Inserting in (27) , we arrive at 

E (O) > (1 - a(^))|i;|ef - 2\n\,a{i) - 2K\Bai\. 

Clearly 

n<M<\l',\ + \Bat\. (29) 
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Using then that \ef^\ < m£Lx(«;, k'), we obtain that 

for some uniform constant C > 1. 

Finally, we use that O, G TZfj to estimate \Bai\. It is clear that there is a 
constant a' > 1 (depending only on A) such that 

Bae C {x G I d{x, di^) < a'i} , 

thus, for Ki := a'C and < d (recall fj is defined on [0, c')), 

\Bai\<mn{T^^<wn 



which ends the proof of Lemma 5. □ 
Step 3. Convergence of ef and inffje^. Let us introduce 

where we recall that ee{g) = E{geA)\iA\-'^ . Notice that by (A3) (and 
taking L — > oo), 

av \ m 

By Lemma 5 (applied with fj = rj and Q = i'A), we have for any £,£'>!, 

ef - ai(£) - Kir/ {k^IIH!) < ef - ai{£) - mrj (k2^/^') < < e^7 (30) 

where K2 = |A|~-^/^ and provided K,iK2i/i' < c. Let (in) and (i'^) be two 
sequences such that 

lim e^ =liminfei? and lim e^T =limsupeiP. 



m . 
av 



Extracting subsequences, we can assume that £n/^n — as n — ^ oo. By (30), 

we deduce that lim supf_^^ < lim inf ^^oo and thus that lim^^oo e 
e exists. By the same argument and (30), it coincides with lim^-^ooej 
Therefore, we have proved that there exists e such that 

lim ef = lim ef = e. (31) 
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The same argument applied to 



:= inf eiMA) < 
ReSO{3) 



shows 



Hm ef = e. (32) 

»oo 

We end this step by proving that ei{R) := e{RiA) converges to e in 
L^{S0{3),dR). Indeed, let us introduce gi := — e and write gt = — gj 
with gj,gj > and gjg^ = 0. Then, since gj = g^ x le^<e(.( )<e) we infer 

< 5," < - e| 
which proves that g^ converges uniformly to on SO{3). Now, 

/ 9i= [ 9i + (ef - e) 
JSO{Z) J so (3) 

which converges to as ^ — oo and shows 

lim||e^-e||^i(50(3),(iil) = 0- 

Step 4. Uniform convergence of towards e. Let us fix some (small 
enough) constant 5 G (0, 1). Since A is open and convex, it is clear that 
there exists a neighborhood A of Ox Id in G such that Ug^Agi'f — S)A C A. 
We can choose A of the form A = B{0, r) x W where is a neighborhood 
of Id in SO{3) and r > 0. Also one has that for £q large enough, £qA and 
45(1-5) A both belong to 7^ by (PI) and (P3). Thus for any {g,g') e AxG 
and any £ > Iq, g'iA and g'ig{l — S)A both belong to TZ. Using (P2) and 
(A4), we infer that 

E{g'eA) < E{g'gi{l-S)A) + K\{g'eA)\{g'gi{l-S)A)\ + \eA\ai\eA\), (33) 

for any g' e G, g e Ae := 5(0, r^) x and ^ > £o. Denoting g' = {u', R') 
and integrating g over A£ in (33), we obtain 



JB{u',re)xR'w FlJ- - <^j^l 

By (A3) and Lemma 4 

n^r^f dXig fun'-Jl'^^ -f e{m{l-5)A)dR 

< \R'W\a{r£), (35) 
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and 

\R'W\-^ [ dRe{m{l-S)A)-e 

J R'W 

(36) 

which converges to when £ — ^ oo as proved in the previous step (here 
denotes the Haar measure of the set R'W in SO {3)). Thus 

supe^ <e + a{r£) + \W\-^ \\ei{i-6) - e||ii(50(3),dR) + + "d^^l)' 

Since we have aheady proved that Hm^^oo inffj = e in Step 3, Theorem 1 
is proved by first passing to the Umit as £ — oo and then S 0. 

Remark 15. Adapting the previous arguments, one proves that e^(-R) = 
e{R£A) converges to e uniformly on 50(3). 



< \W\ ^\\ei{i-S)-e\\mso{3),dR) 



Step 5. The limit e does not depend on the reference set A. Assume that 
there exists another convex set A' with the same properties (PI) (P3), 
that E also satisfies (Al) (A5) with A', and that M5 = Ma U A^a'- 
Denote by e' the limit of E{LA')\LA'\-^ . Applying (24), we obtain 



ILA'I - ' ' "VL|A'|1/3 

hence taking two sequences L„ ^ 00 and in 00 such that in/ Ln ^ 0, we 
see that e' > e. The other inequality is obtained by interchanging A and 
A'. 



2.2 Proof of Theorem 2: limit for general domains 

Let us fix a sequence {On} ^ "^j such that — 00 and diam(On)|0„|'~^/^ < 
C. 

Step 1. Lower hound. We start by proving that 

liminf:^^^>e. (37) 

This is indeed a simple application of Lemma 5 with $7 = J7„ G C = 
and £ := |J^„|^/^. Applying (24), we obtain 
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which proves (37) since ej^^^ |i/6 — e by the proof of Theorem 1. The rest of 
the proof is then devoted to showing the upper bound 

limsup:^^<e. (38) 

Step 2. The inner approximation ofQn- By assumption, we have that the 
smallest ball Bn containing 0„ satisfies > 6' = independently 

of n. On the other hand, since A is open, it contains a ball B{0,r), r > 0. 
Translating and dilating A, we therefore obtain that there exists A[j := 
dnQn^ for some {gn\ C G and — oo such that i7„ C A^ for all n > 1 and 

l>i^>c>0. 

I nl 

We now consider the tiling {/xA}juer of M^. We fix a sequence {^n}n>i 
satisfying in/^'n ~^ ^ ^'^^ which will be specified below. We apply (A6) 
with Q, = A^ which belongs to TZ by (P3). 

We notice that there exists a G-invariant (with the obvious G action on 
G/V) measure d\ on G/V such that 

/ f{g)dX{g) = [ />)dA(a), (39) 
J Jg/t 

where to any function / € L^{G) we have defined / : G/F ^ C by 

/([5]) = E/(5/^)' 

g being any representative for the left coset [g] = {gfi, /x G F} G G/F. Taking 
f{g) = t^igx) for some fixed x G and using that A defines a F-tiling, 
we obtain in particular from (8) that 

[ dX{a) = |A|. (40) 
Jg/t 

Hence (A6.3) can be equivalently written 

^ f dX{[g]) E lt}i9f^'9i^) ^ -iKH^n). 
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We deduce from (40) that there exists a rotation and translation G G of 
the tihng such that 

J2 lt:i9ntJi,9n^) > -iKH^n) > -\n^\c-^a{e^). (41) 

Let us introduce the notation A,i(/x) := ingnfj{^ + T{in))^- We define the 
set of /Lt's which are such that the associated A„(/x)'s are strictly inside 

T'n := {/w e r I An{^^) C Qn and d{An{^^) , dQn) > 6} 

where 5 is the constant appearing in (A4). We finally introduce the inner 
approximation of ri„ by the tiling (see (18)) 

An = A^(^„),^„,p„(J^n) = U C 

By Proposition 2 and (P4), we know that An & TZ' C. TZ^i for in large 
enough (notice that we use £n/^n ~^ 0)- We can thus apply (A4) and obtain 

Ei^n) < E{An) + K\nn \ A^] + |0„|a(|0„|) = E{An) + o(|r2„|) (42) 

and it remains to estimate E{An). 

Step 3. Estimate on E{An). For the sake of clarity, we introduce the 
notations 

A' A' A' 

En (/^) = {gnlJ-) , In{lJ;l^) = I^J" {gnl^, Qnl') and Sn (V) = S^J {gn , V) . 

By (A6.2), we have 

An \oi((.rn 

)• (43) 

Then, by (A6.1), we infer 

E{An) < E{A'^) - Enitx)-^ Inif^,l^)+Sn(r)-SniPn) 

/xer\-p„ //er,t/er\p„ 
-\ H /n(/i,z^) + 2|A;|a(£J 



<e(a;)- E ^"(/^'^) 

I Y Inil^,i^) + 2\K\aiin) (44) 



2 
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Figure 4: The three sets An Q 0,n Q and the tihng {A„(/u)}^ used in 
the proof (with r = here). 

where we have used that Sn(r) — Sn('Pn) < by (21). 

Notice by Theorem 1, E{A'^) = e| A^| + o(| A^|). In order to estimate 
the terms of the right hand side of (44), we shah need the following two 
lemmas, whose proofs will be postponed until the end of the proof of this 
section. 



Lemma 6. There exists a function a2 tending to at infinity, such that for 
any P C 



Lemma 7. There exists a decreasing function 03 tending to at infinity, 
such that for any V, V C T, 



We now show how to end the proof of Theorem 2 assuming that Lemmas 



{V)\+ {Eni^i) - e| a; n A„(/i)|) < |A;|a2(4). (45) 




(46) 
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6 and 7 hold. Indeed we deduce from (44) and Theorem 1 that 

E{An) < e\K\- J2 e|A;nA4/x)| + |A;|«3(4)ff^y 

+|A;|(a2(4)+2a(4)) + o(|A;|) 

< e|J^„| + |A;|a3(£„) (1^) 

+\K\iM^n) + 2a{£n)) + o{\A'J). (47) 

Now, we can find a sequence {in} satisfying in oo, in/i'n ~^ 
«2(^n)(-^n/^n)^ — > as n ^ OO. It Suffices to take for in the solution of the 
equation 



which always exists and satisfies the desired properties since 03 is decreasing 
and tends to at infinity. Then, the upper bound (38) is a consequence of 
(42) and (47) and it only remains to prove Lemmas 6 and 7. 

Step 4. Proof of Lemmas 6 and 7. We shall need a special treatment for 
the An{fi) which are close to the boundary of A^. For this reason, we first 
introduce the following set 

Bn := {^^ G r, I n A„(^)| < \An{fi)\C{in)] 

where ( is some function tending to at infinity. Indeed, the following choice 
will be convenient: 



din) = max (^ai (|A|V6^y^' , |A|V6^-V2^ . 



(48) 



We start by proving the 



Lemma 8. There exists a function 04 tending to at infinity, such that the 
following holds, for n large enough: 

w cr\« i^(A; n (Anifi) U A„(^.))) ^ _ 

^^^eT\Bn,.eT, |a; n (A„(;.) u A„(.))| ^^-"^(^")- (49) 

Proof of Lemma 8. We first notice that when G F \ B„ and G F, then 

|A; n {Aniix) U An{v))\ > \An{ix)mn) (50) 
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by the definition of B„. We now prove that n (A„(/i) U A„(i/)) satisfies 
the r/n-boundary property for the function r]n{t) = cC,{(^n)~'^^^t where c is a 
constant. For the sake of simphcity, we only treat the case where fj, = v E 
r \ Bn, the argument being the same for the general case. 

Since A is a polyhedron, then 9 := A^ n A„(/x) is also a polyhedron. 
Therefore, for any a > 0, there exists a uniform constant </ such that 

t/i^ < a ^ |{x I d(x, de) <t}\< c'{ln)H. (51) 

To prove this, it suffices to notice that the boundary of Q is composed of a 
number of faces which is at most twice the number of faces of A, and that 
the area of each of these faces is at most of order {InY- If we assume that 
< 1, then since |G| < |A„(^)|, we have t/in ^ o for some a and we 
can apply (51). We get that when t|9|~^/^ < 1, 

where we have used that |0| > | A„(/x)|C(^n) since jj, ^ Bn- This proves that 
satisfies the ry^-regular property with rin{t) = c({in)~'^^^t, t G [0, 1). 

We are now able to apply Lemma 5 with $7 = C A„(^) and I = 
|e|V6|A|V6. Note that = < ln^^^C{in)~^^^ which 

tends to zero as n — > oo by the definition of C(^ra)- Hence for n large enough 
we have < l/^i and we can apply Lemma 5. Denoting /3(r) = e—ef 

which tends to when r — oo by the proof of Theorem 1, we obtain 

^ > e--/3(|e|V6|A|V6)-^^ai(|G|V6|A|V6) 

C(4)-^/='|A|V6 



-KlC- 



> e-/3(|e|V6|A|V6) 



e|i/6 

ai (|A|V6^) c|A|i/6 



C(4) C(4)5/6(^n)V2' 

Using the definition (48) of Ci^n), we obtain 

>g_ " ( „. i\A\7;./fj \ ri(o /iai-\ v. 



sup {/?(*)}-«! (lAI^V^)' -C(V|A|3)" 

III \M1^III M /fi V / 
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l®l t>|A|l/6(^„)l/2^(£„)l/6 

which ends the proof of Lemma 8. □ 
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Proof of Lemma 6. By Theorem 1, we have -E'(A^) = e| A^| + | A^|q;(^^) for 
some function a. tending to zero at infinity. By (A6.1) and the choice of Qn 
to ensure (41), we have 

Enilj) - Sn{T) + \l^'nWn) < E{A'J < e\K\ + I A^l sup a(t) 

<e^|A;nA„(/x)| + |A;|a5(4), (52) 

for some function tending to zero at infinity. In the first Une of (52) we 

have used that in/^'n — 1 (indeed we shall choose in/^'n ~^ 0). In the second 
line, we have used that A is a polyhedron, hence belongs to TZtja some 
VAit) = a\t\, which easily shows that 



^|A;nA„(M)|-|A;| 



< c\A'jT{en). 



(53) 



(54) 



By the nonpositivity (22) of Sn, we deduce from (52) that 

J2 [Eniti) - e|A; n A„(/x)|) < |A;|(«5(4) - a(4)). 

By (A6.2) and (A6.5), we have for all /x G T 



Let us then consider some P C F. Summing over fj, E V and using Lemma 
8 with 11 = u and (A2), we obtain 

Y^^M ^ J2 (e-«4(4))|A;nA„(/x)|-K Yl |A;nA„(/x)| 

> e^|A;nA4M)|-('« + e) Y \K^^n{l^)\ 
-2(a(4) + a4(^n))|A;|, 



where the constant 2 comes from (53). Then, we notice that, for n large 
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enough and by the definition of H„ 



A„(M)n5A;^0 
\ 2 



All this shows that 



^ - e| a; n Aniii)\) > -\A'Jae{£n) 

for some function as independent of V and tending to at infinity. 
Applying this result to 7^ = F and recalling (54) , we deduce that 



(55) 



J2{En{fi)-e\A'^nAn{fi)\) 



(56) 



with 



"7(4) = |a5(4) - a(4)| + |a6(4)|- 
We can now deduce the upper bound for any V by 

^ {Enifi)-e\KnAM\) 

<|A;|a7(4)- (^n(/x)-e|A;nA„(/x)|) 

< |A;,|(a7(4)+a6(4)). 
Also (55) and (52) give 

-s„(r) < |A;|(a5(^n) - a(4) + «6(4))- 
Hence by (22) and (21), 

< -Sn{V) < -Sn{T) < |A;|(a5(4) - a{in) + ae{in)) 
for any V QT. This ends the proof of Lemma 6. 



□ 



Before turning to the proof of Lemma 7, we quote the following general 
property of sf: 
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Lemma 9. Assume that sf : G X {V C T} ^ R satisfies (A6.5) and 
( A6.6) . Then we have for any 7^ C F with < oo 

'ii9,r)<^ s^ig,{n,u}). (57) 

We shall give the proof of Lemma 9 later on and rather turn to the 

Proof of Lemma 7. Let V,!^' C T. By (A6.4), we have /„(//, z^) = if 
An(/i) n A'j = or A„(z/) n A'^ = 0. Hence we can prove Lemma 7 assuming 
that Anifi) n A^ / for all /Lt G 7^ U 7^'. In particular 



3 



#{VuV')<C{f) (58) 



for some constant C. Then we write 



\ii^v jX^V fl^V 



/iePnB„,i/6P'ne„ 



By (A6.2), we have for any G L, 



z^) - Sn({/^, z^}) > ^(A; n (A„(/x) U A„(z.))) - £;„(/x) - En{u) 

- |A; n (A„(m) U A„(z/))|a(£„). (60) 

If /X, G Bn, we have by (A2) 

> -(K + a(4))|A;,n(A„(Ai)UA„(z.))| -^„(/i)-^„(z>) 

> (e|A; n A„(//)| - En{^i)) + (e|A; n A„(z/)| - En{v)) - C(4)'C(^n)- 

(61) 
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Tthus by Lemma 6 



> 



ij.evnB„, veV'nBn 



i' 



(62) 

where we have introduced 
and used (A6.4). 

If now n£V\Bnandu e V, for some V,V' CT, we use (A6.2) 

and Lemma 8 to obtain 

In{l^,v) - sn{{li,y}) > ejA^n (A„(/x)U A„(z/))| 

-|A; n (A„(/x) U A„(z/))|(a(4) + 04(4)) 
> (e|A; n A„(/x)| - EM) + (e| a; n A„(/x)| - E^) 
-|A; n (A„(/x) U An{v))\{a{ln) + ^4(4)) 
-as (4) I A„ (/x) I (5n (/x) (5„ (z/) 

for some function ag tending to at infinity and wlicre we have introduced 

r 1 if A„(^)nA;7^0 

"^^^ ~ \ otherwise. 

Summing then over /j, E V \Bn and v E V' and using (A6.4) and Lemma 
6, we obtain 



a9(£„). (63) 



n&v\Bu,vev' iJ.ev\Bn,veP' 
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By the nonpositivity of s„ and Lemma 9, we have 



ix^v \x^v 

>3 sn{{^iM)>mv^v')sn{vv^r')>-c(^f] |A;|a2(4). 



(64) 



where we have used (58) and Lemma 6 in the last estimate. We conclude 
the proof of Lemma 7 by (59), (62), (63) and (64). □ 

We now give the 

Proof of Lemma 9. For simplicity we use the shorthand notation s{V) = 
s^{g,V). We prove (57) by induction on #r. If #V = 2, then (57) is an 
identity. Let us assume that (57) holds triic for any V CV with < N 
and prove it for some V with = + 1. 

By the strong subadditivity (A6.6), we have for all pairs {lJ.,i^} C V 
with fJ- ^ v 

s{r) < s{v \ M) + s{{fi, i^})<^ Y ^({"' /^}) + ^({/^' ^})- 

a,per\{ij.} 

Now we sum over all pairs {/x, u} C V: 

^^^^^^ s{r) <^Y E ^({«'/^}) + ^ E ^({^'^})- (65) 

Notice 

Y Y ^({«'/^}) = E E E 

= N{N-l) Y «({«,/?})• 

a,l3eV 

Inserting in (65) we obtain 



35 



which proves (57) for #7^ = iV + 1. □ 

2.3 Proof of Proposition 2: regulcirity of inner approxima- 
tion 

Let f2 be a fixed set in TZj^, g e G, t e [0, tq] and define 

A= (J A(/x) with A{n) =egti{l + T)A. (66) 

/isr I A(At)cn, 

d(9A(^t),9n)>5 

We also introduce V := {^J, e T \ A(//) C A, A(/i) n 9^ / 0}. We 
have to estimate |{x | d{x,dA) < \A\^/^t]\. Let us assume first that t < 
|A(/x)|^/^|^|~^/^c. In this case, we write 

{x I d{x,dA) < \A\^/H^ C U |x I d{x,dA{n)) < \A\^/H^ 
and, using that A(/x) G TZr^, we infer 

|{x I d{x,dA) < \A\^/h}\ < (#P)|A(/x)| X r, (^J^^^ (67) 

(recall that 77 is by assumption defined on [0, c)). If we assume that i is large 
enough compared to 6, £ > ig, then there exists a constant 7 such that 

V/x e T', A(/x) C |x I d(a;, dQ) < 7I A(/x)|^/^} . 

In each A(//), we can choose a ball -B(yu) of volume proportional to that of 
A(/x) such that all the -B(/x), fi eV never overlap. This implies 

(#P)|A(/x)| < 7' |{x I d(x,9J^) < 7|A(/x)|V3|| 
for some constant 7'. Using now that Q has an ?7-regular boundary, we have 

\{x I d{x,dn) < 7|A(/x)|V3}| < X ry (^ ^'^ff' ^ (68) 

when 7I A(/x)|^/^|ri|"^/^ < c i.e. when | A(/Lt)|^/^|ri|~^/^ is smaU enough, 
£ < 6o\n\^/^. Hence 



{#V)\A{fi)\<j'\n\xr, 



' 7|A(m)|V3 ' 

|0|V3 
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Inserting in (67) , we deduce that 

|{.|a(.,..,.M|V3.}|,y,„,,(JiJ|L),(«) (-) 

when t < |A(//)|i/3|^|-i/3g^ gy 

same arguments as above, it can be 
proved that \ < \A\ +C\^l\rj (7I A(;[/)|^/^|f]|~^/^) , i.e. we have for instance 
< 2\A\. We now use the specific form of 77 G 'r){u) = a\u\'' to infer 

||x I d{x,dA) < < 2j'a-f''\A\ri{t). (70) 

If on the other hand t > c| A(/Lt)|^/^|74|~^/^, then we just say directly that 
I d(x,5A) <t|y4|V3| ^ |x I d(x,5J^) <7|A(/i)|V3 + f|A|V3| 

C |x I d(x, dn) < (c + 7)t|A|^/2} 

and thus for t < c{c + 

||x I d{x,dA) <t\A\^/^j^ < 2\A\{1 + -ffj] (t) . (71) 

As a conclusion, A e TZfi for fj{t) = mqit) on [0, c/m) and with m = 
max (2(1 + 7)^27'a7^c + 7). □ 



2.4 Proof of Proposition 3: (A6) implies (A5) 

For any > we denote by £ = /(/) the smallest positive number such 
that i' = ^(1 + r(£)). Notice that lim^/^00 /(^') = 00 since £(1 + ||r||^oo) > i' ■ 
For any fixed we apply (A6) with the corresponding ^ = f{l'). 

Note that the right hand side of (19) indeed only depends of [g] G G/V 
and not g G G. Hence we can integrate (A6.1) over G/T and, using both 
(40), (A6.2) and (A6.3), we obtain 



E{n) > ^ / dXi[g]) E{n n ^5/^(1 + r(£))A) - 3\nia{^) 

-a{i)-^ [ dA(b])^|J7n£5Ml + T(^))A|. 
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We then use (39) to infer 

^ / dA(b]) J] |0n^5M(l+rW)A| = ^ / \nne9{i+T{e))A\dX{g) 
I I JG/r I I JG 

= Jl^JJ^^ 5^(1 + r{e))A\dX{g) = \n\{l + r{£)f 
by (17) and similarly 




>j^J^E{i}n ge'A)d\{g) - ac((1 + t(£))'^ - l) 
where we have used (A2). This easily proves (A5). □ 

2.5 Proof of Lemma 1: regulcirity of open convex 

Let us fix some open and bounded convex set K, containing 

\{x G I d{x, dK) < t}\ = \K + tB\ -\K ^ tB\ 

where B is the closed unit ball and where we recall that A ^ B 
b), see [25, Chap. 3]. As K is open and contains 0, there exists 
that rB C K. Hence we have for all < t < r 

K + tB C {l + t/r)K and K tB DK {t/r)K = {1 - t/r)K. 

Therefore 

\{xeR^ \d{x,dn)<t}\ < \{l + t/r)K\-\{l-t/r)K\ 

Of 

= {{1 + t/rf - {l-t/rf)\K\<j\K\. 

□ 
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